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Abstract 

The effect on spectral statistics and on the revival probability of intruder 
states in a random background is analysed numerically and with perturba- 
tive methods. For random coupling the intruder does not affect the GOE 
spectral statistics of the background significantly, while a constant coupling 
causes very strong correlations at short range with a fourth power depen- 
dence of the spectral two-point function at the origin. The revival probability 
is significantly depressed for constant coupling as compared to random cou- 
pling. 
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I. Introduction 

About a quarter of a century ago the effect of an intruder state on the 
spectral statistics of a Gaussian orthogonal ensemble (GOE) was discussed 
[1] in the context of isobaric analogue states in nuclei. The limited numerical 
possibilities revealed a very small effect in short range correlations. Recently 
the problem of intruder states has regained interest in the context of the 
analysis of revival times [2] and of state sensitivity [3] of quantum chaotic 
systems. The original model [1] proposes a random coupling of the intruder 
with the background, quite similar to the one proposed in [3] . Ref . [2] , on the 
other hand, proposes a constant coupling and finds an anomalous behaviour 
in the revival probability. 

In the present paper we shall discuss differences between the two models 
in terms of spectral statistics. The random coupling case is revisited and with 
better statistics the deviations from the statistics of the background seems 
even smaller. Constant coupling on the other hand introduces a quartic re- 
pulsion between levels, which leads for the particular case of the intermediate 
coupling chosen in [2] to a nearest-neighbour spacing distribution that is quite 
similar to the one of a Gaussian simplectic ensemble (GSE); at intermediate 
ranges we find a stiffness of spectral fluctuations that exceeds even the one 
of a GSE. At long ranges the entire effect disappears and we revert to the 
GOE behaviour. 

Next we analyze the revival probability. On average it is significantly 
larger for random coupling than for constant coupling. This is readily under- 
stood in terms of the size of the intensity fluctuations. Furthermore, for the 
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constant coupling case the revival probability shows a peak, that is related 
to the anomalous spectral statistics mentioned above. 



II. The models 

The basic question we wish to address refers to the effect of an intruder 
state immersed into a "chaotic" background of states. This background, 
according to a common assumption, can be modeled by a GOE for a time 
reversal invariant system. The open question refers to the type of coupling 
between the background and the intruder. In the early paper [1] the natural 
assumption was made, that the coupling matrix elements should be drawn 
randomly from a Gaussian distribution, but that there should be a variable 
parameter to change the strength of the coupling. We shall henceforth refer 
to this case as model I. The model used in [3] is quite similar. 

Model II proposed in [2] is even simpler as it chooses the same coupling 
matrix element between the intruder and all background states in the basis 
where the background Hamiltonian is diagonal. 

We specify our models as follows: Let if_ be a n x n GOE matrix. 
Diagonalyse H and subject the spectrum, centred at zero, to the usual 
unfolding procedure [4] normalising the average density to one. We now add 
the intruder state in the centre of the spectrum by adding a row and a column 
to the matrix. We obtain 



where v = (v±, v 2 --.v n ) gives the coupling between the intruder and the back- 
ground states. The two models differ in the choice of the coupling matrix 




(2.1) 
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elements V{. 

As mentioned above, model I requires the coupling matrix elements to be 
chosen from a Gaussian distribution with width v = \fup. Here v indicates 
the coupling strength. Actually our model does not coincide exactly with 
the one discussed in [1] because there the intruder is placed in the original 
matrix. The diagonalisation of the background will not affect the coupling 
because an orthogonal transformation will leave ( in the limit of large n) an 
ensemble of vectors of Gaussian random variables invariant. The unfolding, 
on the other hand, does affect the matrix, but we expect only edge effects, 
which should vanish at large n. 

Model II simply chooses v = constant. At v — 2.08 we are exactly in the 
case of [2]. 

III. Spectral statistics: numerical analysis 

In order to obtain information about spectral statistics we have to di- 
agonalise H and again unfold the resulting spectra. The calculations were 
performed with ensembles of 1000 matrices of dimension 300 x 300 and their 
reliability checked by several calculations with ensembles of 800 x 800 matri- 
ces. After the diagonalisation of the background states the fifty lowest and 
highest states were omitted to avoid finite size effects. 

This was done for both models and various values of v. At small v we 
are in a perturbative regime, for v — 1 we are in a transition zone which in 
model I coincides exactly with the GOE. 

For v = 2.08 we shall find a spreading width T = 27 measured in terms of 
the average level spacing normalised to one. We show the intensities of the 
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intruder in the eigenstates in Fig l.a for model I and in Fig lb for model II. 
The results are given in both cases for an ensemble of 10 matrices. We see that 
the results fluctuate widely, and more so for model I than for model II. We 
may quantify this by calculating the inverse participation ratio P l = Y^=i If- 
Here Jj = |aj| 2 are the intensities of the intruder in the eigenstates as shown in 
Fig. 1 obtained from the amplitudes that result from the diagonalisation. 
We find for models I and II respectively Pj = 0.036 and P l n = 0.019, in 
agreement with the impression we obtain from the figures. 

To obtain a spreading width we smooth the data by performing an ensem- 
ble average and plotting the total intensities of the intruder in states lying 
in a unit energy interval. The result of this procedure is shown in Fig. 2 for 
Model II. The fit with a Lorentzian shows excellent agreement. For model I 
a similar result is obtained, but not shown. 

At some larger value of T, that depends on the size of the matrix, the 
spreading width will reach the dimension of the matrix and edge effects will 
become important. For very large and very small couplings we shall present 
some perturbation results in the next section, that will essentially explain 
our numerical findings. 

For model I the nearest-neighbour spacing distribution remains very close 
to the one of GOE in accordance with the results reported in [1]. Model II 
on the other hand shows a rather surprising development from a distribution 
apparently rather similar to the one of GOE for v = 1 to one resembling a 
GSE for v = 5 as we can appreciate in Fig 3. Fig. 4 shows the same distribu- 
tion for v = 1 and v = 2.08, but restricted to 10 and 30 spacings respectively 
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near the centre of the spectrum, which keeps us inside the corresponding 
spreading width (see Fig. 2). We now find the GSE type behaviour fully 
developed in both cases. Note that the agreement with the GSE is not very 
good so we must be careful not to jump to conclusions. On the other hand 
the disagreement with a GOE is absolutely convincing. 

We next want to check the long range spectral correlations. Here we 
make use of a well known fact: As we add a single row and column to a 
matrix the eigenvalues of the new matrix are intertwined between the ones 
of the old one, as displayed in Fig. 5. Therefore, in both models the long 
range behaviour must approach that of a GOE. Actually we expect to find 
no significant deviation beyond three average level spacings. 

As a consequence of the above exact result we may limit our attention 
to the intermediate region. For this region a study of the Fourier transform 
is most profitable. Fig. 6 shows such transforms of the unfolded spectra for 
the same values of v and the same energy intervals we used for figure 4. Now 
we restrict our attention to correlations inside the spreading width with the 
same cutoff's as above. Again we see a surprising resemblance with the GSE. 
Indeed we find a very marked peak at the value 1 of the transformed variable 
which we shall as usual call time. For shorter times though the correlation 
hole is even stronger than for a GSE thus confirming our guess that the 
coincidence with the properties of this ensemble is limited. 

The number variance shown in Fig. 7 for the strong coupling v — 5 
displays the expected features. At long range it is parallel to the GOE curve. 
At intermediate range we clearly see the transitional behaviour, which we 
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might have expected. At short range it follows the GSE curve [5] rather well 
though the first dip is lower than for the GSE. 

IV. Spectral statistics: Perturbation theory and the quartic 
level repulsion 

From many studies of spectra in the transition region [6,7] we know that 
the power law behaviour of the spectral two-point function at the origin is 
most sensitive to the characteristics of a perturbation. Thus, a GOE per- 
turbed by the slightest GUE, changes its behaviour immediately [6], though 
the effect will be only at very short range. The same was proven for an 
arbitrary symmetric or hermitian matrix perturbed by a GOE or a GUE, 
respectively [7]. 

To understand the level repulsion in model II, we shall present an argu- 
ment that consists of three steps: First we shall show that triplets ( three 
very close lying states) in the background necessarily lead at least to dou- 
blets (two very close lying states) in the perturbed system. Second we shall 
show that existing doublets are quite unstable and have a small chance of 
surviving in model I and none at all in model II. Finally we shall see that new 
doublets will only appear in model I. Thus in model II only such doublets 
will exist as were generated by a triplet in the background. 

The first point is quite simple to show. It again rests upon the fact, used 
in the previous section, that by adding one row and one column to a matrix 
the new eigenvalues must lie between the ones of the original matrix. If a 
triplet exists in the original spectrum one eigenvalue is trapped between the 
first two members of the triplet and another between the second and third 

7 



members. The two trapped eigenvalues necessarily form a doublet in the new 
system. This doublet may or may not form part of a triplet or a quartet in the 
enlarged system, but that is not important in the present context. Recalling 
the well-known result, that we obtain a GSE by eliminating every other state 
from a GOE [8], we find immediately that the probability for a triplet goes 
as x 4 for small level separation x. This simultaneously constitutes an upper 
bound to the power, as every triplet has to cause a doublet. Note that this 
result is non-perturbative. 

The second point involves a longer perturbation argument. We consider 
an energy doublet E ± , E 2 , coupled weakly to the intruder state. Since we ex- 
pect the doublet to couple mainly with itself and, of course with the intruder, 
we may disregard all other states. Thus we consider the matrix 

/ V\ v 2 \ 

v 1 Ei (4.1) 
\v 2 E 2 ) 

where we set < E 1 < E 2 . The secular equation for the eigenvalues of this 
matrix can be written as 

v 2 

X — E\ X — E 2 

We are interested in the eigenvalues near the E\-E 2 doublet. We therefore 
introduce the variables y = X — E\ and 5 = E 2 — E 1 . In these new variables, 
the last equation can be rewritten as 

y + E 1 = v l + ^-. (4.3) 
y y-S 

Now comes the crucial approximation: Since the energies E\ and E 2 are a 
close doublet and since the perturbation due to the intruder is weak, it is 
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clear that neither the eigenvalue close to E\ nor the one close to E 2 will be 
very far from E\. We can therefore always neglect y with respect to E±, so 
that the last equation reduces to 



which yields an equation of second degree, with two solutions y\ and y 2 - For 
these one readily finds 



This equation is easy to interpret: In a first approximation, both E\ and 
E 2 move independently by an amount vf / E iy due to their coupling to the 
intruder state. If this motion does not bring them any nearer to each other 
than \v\V2lE\ I , nothing more happens and this result, obtained from ordinary 
perturbation theory, remains valid. On the other hand, if the eigenvalues are 
closer to each other than \viV 2 /Ei\, then higher orders of perturbation theory 
cause the eigenvalues to interact and they wind up at a distance larger than 
\viv 2 /Ei\. 

Furthermore no doublets can be formed at distances shorter than \v\V 2 /Ei\ 
It follows that, for constant Vi, no new doublets can be formed at all. For 
distributions of the Vi such as the Gaussian, which allow a finite probability 
for very small values of the v iy it is possible to form new doublets from energy 
levels which were separated before the coupling. In the above example, this 
occurs when v\ ^> v 2 . Then E\ moves a great distance towards E 2 and the 
repulsion does not occur unless the distance is extremely small. 




(4.5) 
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Thus in model I new pairs are created and again destroyed in terms of 
typical avoided crossings and we may reasonably expect to retain the linear 
dependence of the GOE two-point function. In model II, on the other hand, 
such avoided crossings will not occur and the fourth power limit which we 
obtained from the triplet will correspond to the actual distribution. 

Note that the same argument can readily be used for an intruder with 
constant coupling to a Poisson, GUE or GSE background, to derive a power 
law of x 1 , x 7 and x 13 , respectively. Points two and three of the above argu- 
ment transfer immediately and the triplet probabilities for these ensembles 
are dominated at short range by the power laws given above. 

Our argument is perturbative and thus basically will hold for small cou- 
plings v or for distances outside the spreading width. Figure 8 shows the 
level spacing distribution for v = 2.08 for both models at energies taken in 
the wings of the spreading range of the intruder, and we find the expected x 
and x 4 behaviour for models I and II, respectively. Fig. 9 shows the corre- 
sponding results for spacings taken outside the spreading. Note that in the 
latter case the expected transition to GOE like behaviour for model II occurs 
at very small level separations. 

The fact that the figures of the previous section confirm this behaviour 
even inside the spreading width, goes beyond the direct scope of our argu- 
ment, but is in keeping with the experience discussed above, that perturba- 
tive features at short range tend to extend their range as the perturbation 
grows. The deviations from GSE at intermediate and long range are equally 
unsurprising as there is no theoretical reason for true GSE behaviour. 
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In the limit of very strong interaction for model II the above argument 
certainly no longer applies. We actually revert to the GOE except at the 
edges of the spectrum. Essentially the only effect of the intruder is to expel 
one state at each end of the spectrum, while leaving the remainder with an 
unchanged statistic. We may see this by applying perturbation theory in the 
opposite sense. That is we consider the singular matrix 

Ho=(l I) (4.6) 



v 

where is a n x n matrix of zeros. We choose this matrix as the un- 
perturbed Hamiltonian and the background Hamiltonian as a perturbation. 
The eigenvalues of the unperturbed Hamiltonian are easily computed to be 
E± = rky/nv and all other eigenvalues are zero. The eigenfunctions for the 
non-zero eigenvalues are {aj}± = {±l/\/2, l/^/(2n), l/^J~{2n)} and the ctj 
of the remaining eigenfunctions must fulfill the condition Yh=i a « = and 
a = 0. The effect of the perturbation (which is now background) on the 
extreme states ± is well described by first order perturbation theory. The 
remaining levels are degenerate and perturbed by a matrix that is very near 
to GOE. Therefore, we expect to recover the GOE statisics near the centre 
of the spectrum. 

V. The revival probability 

Model II was proposed in the context of a revival analysis [2], in which 
a numerical calculation of the revival probability turned out to show some 
disagreement with a simple theoretical estimate, that used GOE spectral 
statistics. We shall proceed to show that such discrepancies are not resolved 
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if we use the correct spectral statistics. One reason for these discrepancies will 
become appearent in terms of the inverse participation ratio of the intensities. 
We recall that the revival probability for the intruder state i/j(Q) is defined 

as 

p(t) = \(m\m)\ 2 - (5.i) 

If the amplitude of the intruder in the ith eigenstate of the full system is Oj, 
its time evolution is ip{t) = J2j a j exp(iEj t). We thus obtain for the revival 

p (t) = E \aj\ 2 \a k \ 2 exp(i(Ej - E k )t) = 
E l°i| 4 + S \aj\ 2 \ak\ 2 [exp(i(Ej - E k )t) + c.c] (5.2) 

i j<k 

For an intruder with Gaussian distributed amplitudes, we recover GOE 
properties and the ensemble average for this expresion has been evaluated to 
yield 

< P W> = ^2-f (5 ' 3) 
for t > 0. Here n is the dimension of the matrix and b(t) the usual form factor 
[4]. As we have seen that in model I the intruder does not significantly affect 
statistics, we may hope that we can extend the usefulness of expression (5.3) 
to this model. 

Equation (5.3) is not valid for t — as P(0) = 1.0. Indeed the finite 
spreading width T will dominate the short-time behaviour. We shall heuris- 
tically include the effect of the spreading width T by adding a short-time 
term e~ Tt and by estimating the effective number of participating states. 
For this purpose we assume a smooth distribution of intensities Li along the 
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Lorentzian of width T and calculate the inverse participation ratio J2i |^| 2 ~ 
n e ff. The resulting effective number of levels is n e ff = nT in agreement with 
ref . [2] . We thus obtain 

(P(t)) « e- rt + -L[3 - 6(1-)]. (5.4) 

This result clearly contains errors of the order 1/n for short times but we 
may hope that they are not important as compared to the exponential term. 

We performed numerical calculations both for models I and II by diag- 
onalising the corresponding Hamiltonians for each member of an ensemble 
and evaluating P(t) by using eq. (5.2) directly for a range of values of t. 
The ensemble average was simply performed by averaging the results at each 
value of t. The results for model II agrees with those obtained in ref. [2] 
with other techniques. 

We compare in Fig. 10 the theoretical result of eq. (5.4) with the nu- 
merical one for model I and find excellent agreement even for intermediate 
times, for which this is not guaranteed by our argument. 

The theoretical curve of ref. [2] is also included in the picture and we find 
as reported in this reference a very partial agreement with model II. Note 
that the theoretical result of ref. [2] corresponds to eq. (5.4) if we replace 
the 3 by a 1. This would, in our derivation, correspond to an absence of 
fluctuations in Fig. lb. The n-dependence of corrections due to fluctuations 
is discussed in ref. [2] but they are not evaluated. We thus identify two sources 
of disagreement for model II with theory: First, inside the spreading width 
the form factor is not the one of a GOE and second, the inverse participation 
ratio must be greater than the inverse of the effective number of levels, though 
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smaller than for a GOE. 

The similarity of our result with that of [2] may tempt us to use a semi- 
empirical formula of the same type, where both the inverse participation 
ratio and the form factor are taken from the numerical experiment. Such 
a formula turns out to be in disagreement with the numerical result for 
model II. In particular the very low values of revival near t ~ 0.5 cannot 
be explained. This indicates that the factorised structure of the equation is 
inadequate and is therefore a sign of strong correlations between energies and 
amplitudes that go beyond the general lorentzian shape. Note though that 
the Fourier transform we found inside the spreading width gives a qualitative 
understanding of the maximum in revival near t = 2n. 

VI. Conclusions 

The present paper analysed the effect of an intruder state on a chaotic 
background represented by an unfolded spectrum with GOE fluctuations. 
We consider both a random Gaussian and a constant coupling to this back- 
ground. For the random coupling case we retrieve standard GOE behaviour 
for spectral statistics, inverse participation ratio and revival probability of the 
intruder, by taking into account the spreading of the intruder in a straight- 
forward way. For constant coupling on the other hand, we found inside the 
spreading width a surprising GSE-like behaviour at short distances and a 
GOE behaviour at long range. We could give exact proof for the latter and a 
perturbative derivation for the former. As far as the inverse participation ra- 
tio is concerned we find an intermediate value corresponding to non- vanishing 
fluctuations of the intensities as seen in Fig. lb. A theoretical evaluation is 
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missing. 

Concerning the revival probability of the intruder we obtain some qual- 
itative understanding from the Fourier transforms of the spectrum and the 
inverse participation ratio, but the data definitively indicate important cou- 
pling between intensities and energy eigenvalues. A theoretical understand- 
ing is again missing. 

It is interesting to note that the results of sections III and V are unaf- 
fected if we provide the constant coupling with a random sign. A possible 
candidate for an intruder with near constant coupling but random signs could 
be an isobaric analogue state in a nucleus, where the T< states are isolated 
resonances forming the background, while the T> analogue state is the in- 
truder. A separate study will have to show whether the matrix elements are 
sufficiently constant to cause the short range effect discovered in this paper. 
A further question is whether there are sufficient data to see this effect in 
experiment. 
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FIGURE CAPTIONS 

Fig. 1: The intensity spectrum I { of the intruder plotted against the eigenen- 
ergies for an ensemble of ten background matrices with a) Gaussian random 
coupling and b) constant coupling to the background. The coupling strength 
is in both cases v = 2.08 

Fig. 2: Intensities as in the previous figure, but smoothed over an ensemble 
of 1000 background matrices showing the average value for energy bins of 
size 1 in the normalised spectrum. The results for v — 1.0 are shown by x 
and those for v = 2.08 correspond to *. The dotted line gives a Lorentzian 
fit with width T = 6.5 and the dashed line the same with T = 27.1 

Fig. 3: Nearest-neighbour spacing distributions are shown for v — 1.0 (con- 
tinuous histogramme) and v = 5 (dashed histogramme) using an ensemble of 
1000 background matrices. The long and short dashed lines show the results 
for GOE and GSE respectively. 

Fig. 4: Nearest-neighbour spacing distributions are shown for v = 1.0 re- 
stricted to 10 spacings around zero (dotted histogramme) and for v = 2.08 
(continuous histogramme) for thirty spacings around zero, using an ensem- 
ble of 1000 background matrices. The long and short dashed lines show the 
results for GOE and GSE respectively. 

Fig. 5: Graphical solution of the characteristic equation for the Hamiltonian 
(2.1). The locations of the poles correspond to the unperturbed energies and 
the circles to the energies of the full Hamiltonian. 
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Fig. 6: The full line shows the Fourier transform of the spectrum for model 
II. In a) for v = 1.0 restricted to 10 states around zero; in b) for v = 2.07 
and restricted to 30 states around zero. The transforms for GOE and GSE 
are given by long and short dashed lines respectively. 

Fig. 7: The number variance for v = 5.0 for the whole spectrum ( continuous 
line). The results for GOE and GSE are given by long and short dashed lines 
respectively. 

Fig. 8: The histogrammes give the nearest-neighbour spacing distributions 
for Models I and II taken in the wings of the spreading range. Model I agrees 
well with the GOE result (dashed dotted line) and model II with the GSE 
(dashed double dotted line). 

Fig. 9: Same as Fig. 8 but taken outside the spreading range. Note the 
transition behaviour from GOE to GSE for model II between x — 0.15 and 
x = 0.25. 

Fig. 10: The revival of the intruder calculated according to eq. (5.2) for 
models I and II. The theory according to eq. (5.4) is given by the short 
dashed line and agrees well with model I. The theory of ref. [2] is given as a 
long dashed line. 
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